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Abstract 

SO(IO) supersymmetric grand unified theories [SUSY GUTs] pro- 
vide a beautiful framework for physics beyond the standard model. 
Experimental measurements of the three gauge couplings are con- 
sistent with unification at a scale Mq ~ 3 x 10^^ GeV. In addition 
predictive models for fermion masses and mixing angles have been 
found which fit the low energy data, including the recent data for 
neutrino oscillations. SO(IO) boundary conditions can be tested 
via the spectrum of superparticles. The simplest models also 
predict neutron and proton decay rates. In this paper we discuss 
nucleon decay rates and obtain reasonable upper bounds. A clear 
picture of the allowed SUSY spectra as constrained by nucleon 
decay is presented. 



1 Introduction 



The standard model is unlikely to be a fundamental theory; it contains 19 
arbitrary parameters, 13 of which are the charged fermion masses and mixing 
angles. At least six more parameters are needed to describe neutrino masses 
and mixing (for three active neutrinos). Supersymmetric [SUSY] grand uni- 
fied theories [GUTs] provide a beautiful framework for understanding many 
of the outstanding problems of the standard model^. 

For this framework to be accepted as a description of Nature, the three 
pillars of SUSY GUTs must be verified. These are 

• [I] Gauge and Yukawa Coupling Unification^, |^, ^; 

• [II] Observable Superparticle Spectrum[^; 

• [III] Nucleon Decay[| 

Of these only (I) has been verified. (II) must await Runll at the Fermilab 
Tevatron or the LHC and (III) is in danger of being observed or excluded 
by SuperKamiokande and Soudan II. It is this last feature which is the main 
subject of this letter. In the next section we review the steps in the calculation 
of nucleon decay rates. 

2 Nucleon Decay 

In SUSY GUTs, dimension five baryon and lepton number violating operators 
resulting from the exchange of color triplet Higgsinos dominate; suppressed 
by one power of an effective color triplet mass Mt. Unlike proton decay 
mediated by gauge boson exchange, the value of Mt is a free parameter. It 
is only constrained by requiring perturbative threshold corrections at Mq- 
This constraint, in conjunction with nucleon decay bounds, however, has 
been used to rule out simple SU(5) SUSY GUTs|0]. SO(IO) on the other 
hand, has escaped exclusion, but it is now under siege. 

In equation 1, the Higgs doublets {H^) and the color triplets T (T) 
are contained in a single 10 dimensional representation 10^ of SO(IO). The 
couplings of the color triplet Higgs fields to quarks and leptons are given 
in terms of Yukawa-like 3x3 complex matrices c^g, Cud, Cue, Cqi. These are 
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related by SO(IO) to the Yukawa matrices Yu, Yd, Ye (eqn. 1).[] In a predictive 
SUSY GUT the arbitrary parameters in the Yukawa matrices, defined at the 
GUT scale, are fixed when fitting charged fermion masses and mixing angles 
at low energies. Hence the parameters in the fiavor matrices Cgg, Cud, Cue, Cqi 
are also fixed once charged fermion masses are fit.0 

HuQYjJ + HdQYdD + HdLYeE 
+Q\cgqQT + QcgiLT + UcudDT + UcueET (1) 

Below the GUT scale the color triplets are integrated out of the theory 
giving the dimension five operators (eqn. 2). All dimensionless (dimension- 
ful) parameters are then renormalized to M^.^ We use universal squark and 
slepton masses (mo), gaugino masses (M(i/2)) and non-universal Higgs masses 
{mn^, mnj at M^. 



HuQY^U + HdQYdD + HdLYeE 
jfQ\c,,Q Qc,iL + ^ 



+irQ\cqiQ QcqiL + irUcudD UcueE (2) 



In the effective theory below Mz the coefficients of the effective (dimen- 
sion six) four fermi baryon and lepton number violating operators are deter- 
mined. These are obtained via one loop graphs with squark, slepton and 
gaugino intermediate lines. In general there are LLLL, LLRR and RRRR 
operators generated via gluino and chargino exchanges.^ 

We then renormalize the four fermi operators from Mz to 1 GeV using 
QCD with the multiplicative factor A3lp.OII.pl 



^For more details on notation, see 

^ In our analysis the Yukawa matrices and gauge couplings are renormalized from Mq 
to Mz using two loop SUSY renormalization group equations [RGBs]. We then use a 
global analysis to fit the data[H. This analysis self consistently checks for electroweak 
symmetry breaking and includes the relevant one loop threshold corrections at Mz- 

•^In a more detailed calculation it may be appropriate to have a hierarchy of effective 
field theories to take into account the hierarchy of SUSY particles between 100 GeV and 
3 TeV. 

"^We have not included the contribution of neutralino loops. For more details, see for 
example 1^. 

^Note, A3 is different than Al which appears in ||l^ and is used in many other works 
on proton decay. This is because Al takes into account two different effects, the QCD 
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^ f^i^)^'/'Vi^)^'/''^(ii#4)^'/''^ =1.32 (3) 

The final step is to evaluate tlie matrix elements of these four fermi opera- 
tors between a nucleon and the lepton + meson final state. This requires lat- 
tice gauge theory calculations and usually chiral Lagrangian analysis |r2|, . 
The decay rates depend significantly on the chiral Lagrangian factors aiat 
and Plat where 

PiatU{k) = e,^^<0|Kc/^)M^|proton(k)>, 

aiatUCk) = eaf3'r<0\{u*"d* ^)u'\pToton(k)> 

and U{k.) is the left handed component of the proton's wavefunction. 

We finally obtain the general amplitude for Higgsino mediated nucleon 
decay given schematically by the formula 

Tex A3 (FF) (LF) Mr' Aai (4) 

where 

• FF is a Flavor Factor depending on Yukawa, gauge couplings and 
Cqq Cqi, Cud Cue cvaluatcd at Mz and the specific nucleon decay mode. 

• LF is a Loop Factor depending on gaugino, squark and slepton masses 
roughly as M(i/2)/'"^o ^'^^ ^(1/2) << ^o- 

The Flavor Factor (FF) is model dependent. In our analysis we use a 
particular SO(IO) SUSY GUT with a U{2) x f/(l)" family symmetry with 
the Yukawa matrices given in Jl^. Note, however, that the Yukawa matrices 



running of the dimension 6 operators from Mz to 1 GeV and the running of the quark 
masses from low energies to Mz in order to use the correct Yukawa coupUngs at the weak 
scale. In our analysis A3, including only the QCD running of the dimension 6 operator, is 
the appropriate factor to use, since we are already using the Yukawa couplings evaluated 
at Mz- The general solution to the RGE for the coefficient C(/i) of the dimension six 
operator is given by C(/z) = C{fj,o){-^flj^)^^^''°^ where 60 = H - Infiavors- Note also 



that analytically we have the relation Al = ^3^'^- Given the numerical value we find for 
A3 = 1.32, we then obtain the value of Al = 0.43. For some reason, this differs from the 
value of A]^ — 0.22 used in previous works. We believe this must be a numerical error. 



4 



are fixed to fit the low energy quark and lepton masses and mixing angles and 
the rate for b — > 57. The model is also constrained by electroweak symmetry 
breaking. We estimate, by comparing two quite different models, that the 
model dependence is probably no more than an order of magnitude in the 
rate. 

The largest uncertainties in the proton lifetime, however, enter through 
the value of the effective color triplet mass Mf, through the sparticle spectrum 
via the Loop Factor (LF), and through the strong interaction matrix element 
Plat- We address these uncertainties below. 

Effective Color Triplet Mass - Mj 

The color triplets are required to be heavy with mass of order Mg since 
they contribute to nucleon decay; Higgs doublets on the other hand must 
have mass of order the weak scale. In SO(IO) the Higgs doublets {Hu, Hd) 
and triplets (T, T) are in the field lOjy which is the only Higgs-like rep- 
resentation coupling to standard model fermions. A simple mechanism for 
accomplishing this doublet-triplet splitting exists in SO(IO), known as the 
Dimopoulos-Wilczek [DW] mechanism ||T^. In this mechanism an adjoint (45 
dimensional) scalar obtains a vev of order (B - L) xMg, with B (Baryon) and 
L (Lepton) number and gives mass to the Higgs field IOh- Since color triplet 
Higgs fields have non-zero B-L charge, while Higgs doublets have zero charge, 
only the triplets obtain mass. A simple variation of the original mechanism 
also allows for the possibility of obtaining large (or small) tan/5 ~ 50(2) so- 
lutions with the addition of four fields, ip; ip', i' (16, 16 representations 
of SO (10)) with the unprimed fields obtaining vevs in the "right-handed neu- 
trino" directions. For self-consistency, ip, ip get mass of order Mq and their 



mass and vevs are generated in the SO(IO) breaking sector of the theory[17 
The Higgs doublet and triplet mass matrices are given by 

10^ {^'f 

I (45) W 101, \ 
M(,, ,) = (45) X (^) 10^ (5) 

V (V^) M )\ {i^'f ) 

where the superscript indicates the SU(5) content of the field. We take the 
vevs (45) {B- L)Mg, M Mg and X ~ 10"^ Mg- We then consider two 
cases 
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1. {ip) ^ (7/;) = 

2. (^) ^ {^) ~ O.IMg 

In both cases 1 and 2, the effective color triplet mass (eqn. 2) is given by 

l/Mt = (M-i)ii = XM/detMt ~ X/M^ ~ (10^^ GeV)"^ (6) 

where Mt is the color triplet mass matrix. Note there is actually no color 
triplet with mass greater than Mpianck [0- 
Consider the light Higgs doublets. 

1. Hu, Hd C 10 H are identified as the light Higgs doublets. We have 
A5 = At- = ~ A at Mg and tan (3 ~ 50. A is the universal third 
generation Yukawa coupling given in A 10^^ I63 I63. 

2. The light Higgs doublets are identified as and 7 with 7 = 



XM 



« 1. We then have A;, = A,- ~ 7 A << Xt ~ A. 



The lightest Higgs doublet (besides Hu, Hd) has mass of order IQ-^Mg. 

A limit on the value of M^, the effective color triplet Higgs mass, is ob- 
tained by requiring perturbative threshold corrections to gauge coupling uni- 
fication. At one loop the definition of the GUT scale is somewhat arbitrary. 
A particularly convenient choice is to define Mq as the scale where the two 
gauge couplings, a^, i = 1,2, meet. We define ac = o:i{Mg) = 0^2 (M^) and 
the relative shift in a^lMc) is given by 

63 = (asiMc) - aG)/aG- (7) 

In general, a value of €3 ~ —(2 — 4%) is needed to obtain ~ 0.119. The 
one loop threshold correction coming solely from the Higgs sector (5 + 5) is 
given by||20[ 



e,{mggs)= l^log(^7) (8) 

which is valid in either case 1 or 2 above, if we let 7 = 

Since the Higgs contribution to 63 is always positive, we must therefore 
have a negative contribution coming from the rest of the GUT sector of the 
theory. If we demand that the maximum allowed threshold correction from 
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the GUT sector is —10% (—8%), we then have at most a positive 6% (4%) 
contribution from the Higgs sector (assuming we need es ~ —4%). This gives 
an upper bound on the allowed values of Mt^. 

Note, in the small tan/? regime there is no explicit suppression factor 
entering the coefficient of the dimension 5 operators, since they are all pro- 
portional to y? jMt. The difference in the small vs large tan/? regimes is the 
maximum allowed value of Mt consistent with perturbative threshold correc- 
tions, g We find the bound < 8 x 10^^7-^ GeV (6 x 10^^7-^) GeV for 
e^lHiggs) < 6% (4%). Unfortunately this bound is exponentially sensitive 
to the assumed maximum allowed correction e^{Higgs). 

Finally we warn the reader that, with additional SO(IO) adjoints and a 
clever modification of the DW doublet-triplet sphtting sector, it is indeed 
possible to suppress proton decay via dimension five operators entirely, see 
for example Chacko and Mohapatra[]T7|. In this case, we would unfortunately 
lose a significant test of SUSY GUTs. Moreover such elaborate constructs 
for suppressing nucleon decay seem entirely contrived and unnatural. 



Natural Superparticle Spectrum 

In order for SUSY to provide a solution to the gauge hierarchy problem, 
the SUSY breaking scale Ksusy niust be of order the weak scale. Otherwise 
we must fine tune in order to have Mz ~ rnuiggs « ^susy- In order for 
nucleon decay rates to be consistent with the present data we need to max- 
imize squark and slepton masses, consistent with naturalness, and minimize 
gaugino masses, consistent with present experimental bounds. 

How heavy can squarks and sleptons be? Since the first two families of 
squarks and sleptons couple weakly to the Higgs bosons, it has been argued 
that it is still natural to have heavy first and second generation squarks 
and sleptons as long as the third generation squarks and sleptons are lighter 
than, say, one TeV[^. In fact in SUSY SO(IO) with Yukawa unification and 
SO(IO) boundary conditions at Mq for soft SUSY breaking mass parameters, 
it was noted that the third generation squarks and sleptons are naturally 
lighter than the first two generations due to RGE running p2|, ^. This is 



due to the fact that in this limit, Af, = Ar = A* ~ 1 has the effect of driving 



"This result was also discussed in [[L9 
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the third generation scalars to lower masses. For mo = 3000 GeV and large 
tan/3 we find all third generation squarks and sleptons are lighter than 1 
TeV, except for the left-handed stau and tau sneutrino. These have mass 
~ 2 TeV. If we estimate the contribution to 5m|^ ~ ^^(m^ + rhl^) and 
demand 6m'jj < (130 GeV)^, we find ^ m^^ < 2100 GeV. Thus we avoid 
fine-tuning in the effective theory at the weak scale. 

Note however that for mo = 3000 GeV we must still do some fine-tuning in 
order to obtain the correct electroweak symmetry breaking via RGE running 
from Mgut to Mz- This is because all three families of squarks and sleptons 
have large mass during most of the running. This fine tuning can be avoided, 
however, if we take large masses at Mq for the first two families of squarks 
and sleptons, while keeping the mass of the third family less than ~ 1 TeV. 
This is certainly consistent with the SO(IO) GUT x U (2) family symmetry 
which we are considering. It may even be possible in this case to increase 
the squark and slepton masses of the first two families above 3 TeV.[] 

For small tan/3, only the stop squark mass is naturally light. Thus, in 
this case, the upper bound on mo is ~ 1000 GeV. 



Lattice Results 



There have been several lattice calculations of the chiral Lagrangian parame- 
ters a^^t, I3lat\ 

-Otlat 



A recent lattice calculation p3l on a significantly larger 
0.015 GeV^. The statistical uncertainties in this 
However systematic uncertainties con- 



lattice gives jSiat 

result are small (±1 in the last digit 
nected with the chiral Lagrangian approach and the quenched approximation 
may be significant, perhaps as large as 50%. Note that the proton lifetime 
is also sensitive to the relative magnitude of aiat and jSiat- In our results we 



use the central value quoted in pSi . 



3 Results 

The process p — > K^v is the dominant decay mode for the proton. In all 
cases we find the rate for n K^v dominates over p K~^u; typically by 

^In order to check these possibihties, we must include two loop RGE running as em- 
phasized by Arkani-Hamed and Murayama . We will investigate this further elsewhere. 
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a factor of 2 - 4. However the best experimental bound is on the latter — 
Super Kamiokande 90% CL bounds on nucleon decay based on 61-ktonyear 
exposure |T^ 



t{p K+u) > 1.9 X 10^3 years (9) 

and we use it to establish that SUSY GUTs are under siege. 

The result of this analysis is the theoretical upper bound on the proton 
lifetime given by 

T{p -> K-'u) = 4.7 X 10^^ years x C-^^l^f x { .^J/^^y f (10) 
for mo = 3000 GeV, M(i/2) = 175 GeV and tan/3 ~ 54 or 

T{p K^u) = 1.0 X 10^^ years x C-^^^l^f x {^^Jh^f (11) 

for mo = 1000 GeV, M(i/2) = 300 GeV and tan/3 = 2.0 

Let us briefly summarize the theoretical input to these upper bounds. 

• This upper bound assumes a very conservative upper limit on the color 
triplet Higgsino mass Mf, obtained by demanding that the Higgs con- 
tribution to the one loop threshold correction to gauge coupling unifi- 
cation at Mc is at most 6%. 

• Gaugino masses are taken to be near the allowed experimental lower 
bounds. While squark and slepton masses must be near the upper 
bounds allowed by naturalness. In SO(IO) with universal squark and 
slepton mass mo, we have mo ~ 3000 GeV for tan /? ~ 50 or mo ~ 1000 
GeV for tan/3 ~ 2. 

• We use the central value for jSiat given in |2^. If the reader prefers 



another value, the new bound can be obtained by a simple rescaling. 

• Any further uncertainty depends on the specific model for the Higgs 
doublet and triplet Yukawa couplings; we estimate this uncertainty to 
be at most an order of magnitude in the lifetime. The upper bound 
given here is for the particular SO(IO) SUSY GUT with a U{2) x f/(l)" 
family symmetry which provides an excellent fit to charged fermion 
masses and mixing angles |]T^ . 

^FoT esiHiggs) < 4% simply change the value of Mt in eqn. |l^ (|ll|) to 6x 10i^(4x lO^O). 
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Clearly we have pushed most of the parameters to (or perhaps beyond) 
what the reader may consider reasonable upper or lower bounds. Neverthe- 
less, with these exceptionally conservative bounds we are barely consistent 
with the latest Super-Kamiokande limits on p — K^v (eqn. |^) |jl5| . 



We have shown that the recent SuperKamiokande bounds on proton decay 
severely constrain SO(IO) SUSY GUTs. Recall, simple SU(5) SUSY GUTs 
have already been excluded by this data0]. Some general conclusions may 
be drawn from our analysis. 

• Based on these results the first two generations of squarks and sleptons 
must have mass significantly greater than 1 TeV; while the third gen- 
eration squarks and sleptons can be lighter than 1 TeV, but not lighter 
than ~ 400 GeV. Thus, at best we expect only the third generation 
squarks and sleptons to be visible at LHC. 

• Clearly proton decay must be seen soon IF minimal SUSY GUTs are 
the correct description of nature. 
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